Estimation non-paramétrique de la densité spectrale d'un 
processus gaussien échantillonné aléatoirement 
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Résumé: 

En utilisant une analyse en ondelette, on construit un estimateur non-paramétrique de la densité spec- 
trale d'un processus gaussien à accroissements stationnaires. Dans un premier temps, on considère le cas 



00 . "idéal" de l'observation d'une trajectoire en temps continu et on donne un théorème de la limite central 

T— I ■ ponctuel et une estimation de l'erreur quadratique moyenne intégrée. Puis, afin de mieux correspondre aux 

CN ' applications, on construit un second estimateur à partir de l'observation d'une trajectoire du processus à 
O ' 

(yQ ' des instants discrets, aléatoires et irrégulièrement espacés. Cet estimateur est obtenu en remplaçant les 



coefficients d'ondelette par leurs discrétisations. Nous donnons un second théorème de la limite central, 
avec une vitesse différente. Des simulations et une application à la modélisation du rythme cardiaque des 
marathoniens sont présentées. 

Abstract (english): 

Erom a wavelet analysis, one dérives a nonparametrical estimator for the spectral density of a Gaussian 
process with stationary incréments. Eirst, the idealistic case of a continuons time path of the process is 
considered. A punctual Central Limit Theorem (CLT) and an estimation of the Mean Integrate Square 
Error (MISE) are established. Next, to fit the applications, one considers the case where one observes 
a path at random times. One built a second estimator obtained by replacing the wavelet coefficients by 
their discretizations. A second CLT and the corresponding estimation of the MISE are provided. Einally, 
simulation results and an application on the heartbeat time séries of marathon runners are presented. 

1 Introduction 

Soit X = {X{t),t € M} un processus à temps continu gaussien centré à accroissements stationnaires II 
admet une représentation (voir Yaglom, 1958) 

X(t) = [ (e^*« - 1) • /^/^(O dW(n, pour tout t € M, (1) 
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où W{dx) est une mesure de Wiener complexe adaptée afin que X{t) soit réel et / est une fonction paire, 
positive appelée densité spectrale. On observe une trajectoire (Xt^ , • • • , X^^ ) du processus X à des instants 
ti, i = 0, . . . ,n irrégulièrement espacés ou aléatoires. Le but de cet article est de proposer un estimateur 
de la densité spectrale /(Ç) pour une fréquence < Ç < cxd et d'étudier sa convergence. 

A titre d'exemple paradigmatique, si /(Ç) = |Ç|~(2^^+i) pour tout Ç G M* , alors le processus X est 
un mouvement brownien fractionnaire (mbf) d'indice de Hurst H et de paramètre d'échelle a; la fonction 
de densité spectrale est entièrement déterminée par les deux seuls paramètres (H, a). On se trouve dans 
un cadre d'estimation paramétrique de la densité spectrale /. De plus, la fonction ln/(Ç) est une fonction 
affine de ln(|Ç|), ce qui correspond à l'autosimilarité du processus. 

L' autosimilarité est très souvent supposée sans vraiment être statistiquement vérifiée ou théoriquement 
justifiée. Dans cet article, nous nous intéressons à l'estimation de la fonction /(.Ç) dans un cadre non- 
paramétrique. Cette approche peut permettre ensuite de vérifier graphiquement l'autosimilarité, ou, par 
exemple, de déterminer une ou plusieurs bandes de fréquences de fréquences ]wi,u;î+i[ telle que /(Ç) = 

^2 . |Ç|-(2H.+1) pQ^^ a;^ < ^ < uj,+i. 

Le plan de cet article est le suivant. Dans la section 2, nous précisons le modèle statistique et les 
hypothèses utilisées. Dans la section 3, on montre un théorème de la limite central (TLC) pour l'estimation 
de la fonction de densité spectrale à l'aide d'une analyse en ondelette théorique, c'est-à-dire avec des 
coefficients d'ondelette calculés par une intégrale sur M. Dans la section suivante, nous montrons que 
la discrétisation des coefficients d'ondelette découlant de l'observation à des instants discrets aléatoires 
conduit à établir un second TLC. Dans la dernière section, nous donnons des exemples numériques. 

2 Hypothèses et notations 

Soit X un processus gaussien centré à accroissements stationnaires défini par l'équation ([T]). On suppose 
que sa densité spectrale / vérifie les conditions (Fl) et (F2) ci-dessous: 

(Fl) La fonction f est borélienne, paire positive et vérifie 

[ (lA|Ç|2)./(Ç)dÇ<oo. (2) 

(F2) Sur IR^ la fonction f est càdlàg et en dehors d'un nombre fini de points < cni < ■ ■ ■ < ojk avec 
K E ]N , sa dérivée f vérifie la condition 

[ (1A|Ç|3).|/'(0| dÇ<oo. (3) 

De plus, on a un contrôle du comportement à haute fréquence: il existe deux constantes Co,Cq > 0, 
une fréquence ujk > G et un réel H > 0, tel que pour \^\ > ljk 

m) < Co|Ç|-(^^+^) et f'{x) < C, 

On dispose de l'observation du processus à des instants discrets sur un intervalle [0,T„], i.e. 

{^X{tQ), X{ti), . . . ,X{tn)) est connue avec to = < ti < • • • < tn = 2n- Les instants d'observations ti sont 

aléatoires et irrégulièrement espacés. On aura besoin que l'hypothèse (S (s)) ci-dessous soit vérifiée pour 
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(S(,s)) // existe une suite de nombres réels strictement positifs {ôn)ne]N telle que (5„ — > et une suite de 

n^oo 

variables aléatoires positives identiquement distribuées et indépendantes {Lk)k<^M telle que 

tk+i - tk := ônLk et ]ELk = l pour e {0, 1, . . . , n - 1}. 

De plus, les v.a. {Lk)keiN sont indépendantes de la tribu Tx engendrée par le processus avec Tx '■= 
(T|X(t), t G iR} et il existe deux constantes < < Mg < oo telles que 

iTT-s < ll-^fclls < Mg pour tout k G IN. 

où pour une v.a. Z et un réel a G (0, co], on note\\Z\\a '■= ""^^^ quand 1E{\Z\"') < oo et 

ll'^lloo = sup^gQ [-^(i^)! P^r convention. 

Un estimateur de la densité spectrale /(Ç) pour une fréquence < Ç < oo est construit à partir de l'analyse 
en ondelette du processus X. Soit une "ondelette mère", et pour n'importe quels échelle et décalage 
(a, b) G 1R*^_ x M, on définit respectivement les coefficients d'ondelette "théoriques" (pour une trajectoire à 
temps continu) et "empiriques" (pour une trajectoire à temps discrets) par 

d^{a,b := a'^^^ j i; (^-^^ X{t) dt (5) 
et e^(a,6) := a-^'^Y.(j i' {^) dt) X (t,) . (6) 

Suivant les différents cas examinés par la suite, l'ondelette '\\) vérifiera des hypothèses issues de la famille 
suivante : 

W(m, g, r) ^/^ : IR^ IR est une fonction vérifiant les trois conditions: 

• pour tout entier n < m, / \t'^ ^l) (t)\ dt < oo, / t'^'4}{t) dt = Q et / \t'^^'^ il) {t)\dt < oo; 

JR Jm Jr 

• il existe une constante > telle que pour tout t £ M 

(i + |t|)^- |V'(t)| < c^. 

• il existe une constante C'^ > telle que pour tout £ M 

{i + \^\y.{\^ic)\ + \^'iO\) < c;. 

Cette famille d'hypothèses est également utilisée dans Bardet et Bertrand (2008). Dans cette présentation, 
nous ferons deux hypothèses supplémentaires qui permettent de simplifier les énoncés et les démonstrations : 

(T) • La durée d'observation Tn est déterministe. 

• La transformée de Fourier de l'ondelette mère, ip a un support compact. 

Remarque 2.1 Les différentes hypothèses imposées sont faibles. 

• L'hypothèse (Fl) est équivalente à l'existence de l'intégrale stochastique ([7]), l'hypothèse (F2) autorise 
à la fois des mbf, des processus à bande finie (voir Rachdi, 2004), processus localement fraction- 
naires utilisés par Kammoun et al. (2007) pour modéliser les battements cardiaques ou des processus 

Rprirnnrl pi finr/ipi /'i>nni ) 
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• L'hypothèse S{s) est satisfaite pour tout s < oo dans le cas où les intervalles d'inter- arrivée suivent 
une loi exponentielle (on parlera alors d'échantillonnage poissonnien) et pour s = oo dans le cas 
déterministe ou le cas intervalles où la loi des intervalles d'inter-arrivée a un support inclus dans 
[a5n,bôn] avec < a < b; la condition ôn — > est naturelle. Seule l'hypothèse d'indépendance 

n—*oo 

des instants d'observations par rapport au processus est restrictive pour les applications en finance, 
mais elle est toujours faite, nous renvoyons sur ce point à Aït-Sahalia et Mykland (2008). A notre 
connaissance, l'hypothèse S (s) avec s < oo est nouvelle et tous les travaux sur les mbf ou les diffusions 
observés à des instants irréguliers ou aléatoires supposent vérifiée S(oo), voir par exemple Begyn 
(2005), Blanke et Vial (2001) ou Hayashi et Yoshida (2005). 

• L'hypothèse W{m,q,r) est classique et peu contraignante, voir la discussion ci-dessous. 

• L'hypothèse (T) est purement technique et peut être supprimée au prix de notations plus lourdes. 
Supposer déterministe la durée d'observation Tn est standard. Cela n'empêche pas les temps d'inter- 
arrivée d'être aléatoires comme par exemple dans le cas d'un processus ponctuel. Dans les appli- 
cations, que ce soit en finance lors de l'observation des cours d'une action ou en biologie lors de 
l'observation du rythme cardiaque d'un individu sain ou malade, on a une durée d'observation T„ 
déterministe (8h pour les cours de bourse, Ih ou 24/i pour les battements de coeur). On peut donner 
dépendant un exemple où est aléatoire : Kammoun et al. (2007) ont étudié le rythme cardiaque 
de coureurs de marathon, T„ correspondant alors à la durée du marathon, qui est indubitablement 
aléatoir^. 



3 Analyse en ondelette d'une trajectoire continue 

Qu'appelle-t-on exactement une "ondelette mère"? Il y a en fait un léger abus, car nous n'avons 
pas besoin que la fonction tp et ses "filles" engendrent K?'{]R). Nous utilisons uniquement des hypothèses 
W(m, q, r) qui sont faibles et signifient que ■0 est une fonction bien localisée en temps et en fréquence 
avec un certain nombre de moments nuls, comme l'est par exemple une ondelette de Lemarié-Meyer. La 
première étape consiste à obtenir une représentation fréquentielle des coefficients: 

Proposition 3.1 Soient X un processus gaussien défini par ([ip avec une densité spectrale f vérifiant 
(Fl) et ip une ondelette vérifiant la condition VF(0,0, 0), alors pour tout {a,b) G ]R*f_ x M, le coefficient 
d'ondelette dxia,b) est bien défini par ^ et admet la représentation fréquentielle 

dx{a, b) = V^ [ 5(aa f^'\i) dW{0. (7) 

JIR 



Démonstration. La démonstration résulte de l'application du Théorème de Pubini stochastique, licite 
grâce aux hypothèses. Elle est donnée en détail, tout comme les autres démonstrations des résultats de cet 
article, dans Bardet et Bertrand (2008). ■ 

On en déduit que, pour une échelle a fixée, le processus [d^{a,b)) ^^j^ est stationnaire, gaussien, centré, 
et a pour variance 

I^{a) := / \i){x)\'^ f{x/a)du. 



TTn mQ-fr-Vi r\ct -J-onnic q lui quccÏ nno rlnréo q1qq^-/-»ivci r^Qv /^nn^T-o un mal-z^Vi r\ct fnr^^-Kall q nno rlnréo à r\cin r^vèc rlol-orminicl-o 
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On montre également que, pour a fixé, le processus [d^{a, n'est pas à longue mémoire (alors que le 

processus X peut l'être). Tout estimateur de la variance de d^{a, •) fournit alors une estimation de I^{a). 

Estimation de la variance des coefficients d'ondelette 

Soit [Ni^)k^]j\] une suite d'entiers telle que Nn — > co. Soit n G W*. Pour une famille de décalages 

n— +00 

61 < 62 < • • • < 1 on estime (a) par la variance empirique 

lN^,^{a) := N-^ ^dlia,bk)- (8) 
k=l 

Comment choisir la famille de décalages {bk)k=i,...,Nn ^ 

Dans le cas d'une observation à des instants déterministes régulièrement espacés, l'usage consiste à pren- 
dre comme famille de décalages les instants d'observation, sans que jamais la question du choix ne soit 
posée. Nous n'avons pas réussi pour l'instant à résoudre la question du choix optimal de la famille de 
décalages. Pour simplifier les calculs et la présentation, nous supposons dans la suite que les décalages bk 
sont régulièrement espacés. On montre dans la section suivante que l'erreur de discrétisation est contrôlée 
si < 61 < • • • < bN„ <Tn — Tn avec p £ (3/4, 1). Ceci conduit au choix suivant: 

(Bl) On se fixe un réel p G (3/4, 1), on note Tn '■= — 2Tn et on définit la famille de décalages par 

k — 1 

h ■■= TP + _ X Tn pour k = l,...,Nn. (9) 

Remarquons que r„ = (ft7V„ — ^i)- La proposition suivante donne une condition suffisante pour que 
l'estimateur I]y^^^{a) satisfasse un TLC. 

Proposition 3.2 Soient X un processus gaussien défini par avec une densité spectrale f vérifiant 
(Fl) et (F2), une ondelette vérifiant la condition PF(1, 0,1/2) et (feî)i<i<Ar„ vérifiant^. Alors si 

Tn > OO, Nn > OO et T„/7V„ > 0, 

ra— »oo n— >oo n^oo 

^{lN^,^{a)-Ma)) M(o,iTTa^ [ \i,iax)\^ f{xf dx) . (10) 

Estimation ponctuelle de la densité spectrale 

Dans l'expression de X^(a), la fonction x i— > |V'(x)p joue le rôle de la fenêtre spectrale et on estime T^{a) 
qui est une moyenne de / au voisinage de x = 1/a. Pour estimer / ponctuellement, il suffit de disposer 
d'une suite d'ondelettes V'A vérifiant W(l, 1/2,0) telles que leurs transformées de Fourier convergent 
vers la masse de Dirac en 1. On montre ainsi le lemme suivant: 

Lemme 3.1 Soit ip une ondelette vérifiant la condition W{1, 0,1/2) telle que le support de sa transformée 
de Fourier soit inclus dans un compact [—A, A] avec A > 0. Soit la famille ('0A)AeJR^ telle que pour tous 
A G M*^ et t £ M, ip\{t) = (1/\/A) e** • tp{t/X). De plus, si la fonction f est dans un voisinage du point 
1/a, on a le développement limité quand A — > oo 



/(l/a) + /'(l/a).- + 0((A)-2) 



ni) 
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Démonstration. OnaV'A(Ç) = V^-'0(A(i^ — 1)) pour Ç € M. La première condition de W( 1,0, 1/2) signifie 
que VOv(O) = i^x\0) = et / |t| VaWI dt < oo, or ^(0) = VÂ-V^(-A) = et ^'(0) = A3/2 .^'(-A) = 



pour A > A, les autres conditions se vérifiant aisément. D'autre part 

I^,ia)= [ X\ip{X{x-l))\^f{x/a)dx= [ \^u)\^xf(- + 



Xa 



du, 



puis, en utilisant la formule de Taylor et la compacité du support de ip, on en déduit (jlip . ■ 
Ainsi, pour une suite (A„)„£jv telle que A„ — > oo, on définit l'estimateur ponctuel de la densité spec- 

n—*oo 

traie par /n(Ç) = lNn,ipx„{^~^) /Wi^Wl^- On montre que cet estimateur converge vers /(Ç) avec un biais 
^ ^f'iO 6t une variance équivalente à (Xnt^^) x 7rÇ^^/(Ç) ||^/;||^4/||'(/'|||^2- Si le carré du biais est 
asymptotiquement négligeable devant la variance, on obtient le TLC suivant ainsi qu'un calcul de l'erreur 
quadratique intégrée (MISE): 

Théorème 3.1 Soient X un processus gaussien défini par ([1]) avec une densité spectrale f vérifiant (Fl), 
(F2) et ip une ondelette vérifiant M^(l, 0,1/2) et suppijj C [—A, A]. Pour tout réel a g]1/3, 1[, soit 



fniO 



-2 
L2 



Nn . 

k=l JJR 



i{t - bk 



X{t)dt 



avec la famille de décalages donnée par (0). 



i) Pour toute fréquence ^ g]0, oo[, quand — > oo, Nn — > oo et r^/N^ — > 0, on a 

n^oo ra— >oo n— >oo 



i'-"^/'(7n(Ç)-/(6) ^ AA(0,4vr 



(12) 



ii) Pour toute bande de fréquences finie < < f-^^i < oo, on a le développement limité du MISE: 

4tt 



lE 



fn{î)-fm dC 



UJO 



4 , 
^2 •''^0 



Démonstration. La proposition 13.21 et le lemme 13.11 sont les deux ingrédients principaux pour obtenir 
ce théorème. ■ 



Remarque 3.1 La proposition \3.^ permet d'estimer des paramètres de la densité spectrale dans un cadre 
paramétrique ou semi-paramétrique (par exemple dans le cas du mbf multi- échelle, voir Bardet et Bertrand, 

— 1/2 

2008). La vitesse de convergence est alors enTn , où Tn = bj^n ~bi correspond à l'amplitude de la famille 
de décalages; dans le cadre non-paramétrique du théorème \3.1l la vitesse de convergence est en Tn ^ 
si a = 1/3 + 2e pour un e > 0. Concernant l'estimation du MISE sur une bande finie [a;o,a;i], il n'y a 
aucune difficulté à faire tendre la borne supérieure uJi vers l'infini; en revanche, on a une explosion à basse 
fréauence auand ojn 0. 
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4 Estimation de la densité spectrale pour une trajectoire observée à 
des instants discrets aléatoires 

4.1 Autres méthodes existantes 
Périodogramme empirique 

L'estimation de la densité spectrale d'un processus gaussien stationnaire à partir de l'observation d'une 

trajectoire échantillonnée à des instants irréguliers est un vieux problème en traitement du signal, voir 

Lii & Masry (1994) ou Rachdi (2004) et les références qui y sont citées. Malgré la ressemblance, nous 

ne travaillons pas sur la même question. Tout d'abord, nous considérons des processus à accroissements 

stationnaires au lieu de processus stationnaires, mais ce n'est qu'un détail: dans le cas stationnaire, il 

suffit de remplacer la formule ([I]) par X{t) = / e**^ • f^^'^iC) dW{^). La vraie différence provient de la 

Jm 

modélisation et des applications sous-jacentes: dans les années 1960, il s'agissait de choisir une bonne 
méthode d'échantillonnage afin d'estimer la densité spectrale d'un signal analogique; la perspective est 
aujourd'hui entièrement différente, on dispose de l'observation d'un signal à des instants digitaux aléatoires, 
sans aucune possibilité de les choisir. 

Dans Rachdi (2004) ou Lii & Masry (1994), on choisit un échantillonnage poissonien d'intensité moyenne 
P connue, ainsi que la densité de sa covariance c{u), formules (2.1, 2.2 et 2.3) dans Lii & Masry (1994), ces 
quantités sont ensuite utilisées pour définir la fonction 7(n) dont la transformée de Fourier r(A) intervient 
dans la définition de l'estimateur de la densité spectrale (formule 2.14). Une construction différente, bien 
que toujours basée sur la connaissance de la constante (3 et de la fonction c(u) est proposée dans Rachdi 
(2004). Cette approche paraît difficilement transposable à nos données : en admettant que les temps 
d'inter-arrivée soient poissonniens, il faudrait disposer d'une estimation /3 de /3 et surtout d'une bonne 
estimation c(u) de la fonction c(n) . 

Variations quadratiques 

Une autre méthode populaire pour l'estimation du paramètre de Hurst d'un mbf consiste en l'utilisation des 
variations quadratiques généralisées (voir Guyon & Léon, 1989, ou Istas & Lang, 1997). Pour (a, 6) G M^, 

soit Qxio,,b) := X{b + a) — 2X(b) + X{b — a); on estime alors la variation quadratique d'ordre 2 (ordre 

N 

minimal requis pour notre modèle) par VAr(a) := ''^^Qxio-ibk)'^ ■ Le calcul élémentaire 

fc=i 

= 4 / e'^^ ■ sin2«/2) • /^/^(ç) 
donne une formule de représentation analogue à ([7]). On en déduit que le processus (Qx (a, 6))^^^^ est 

/>oo 

stationnaire, gaussien, centré, de variance Vx{o) '■= 64 / sin^(2'u) f{u/a) du, puis que Yxia) converge 

Jo 

vers Vx(a) selon un TLC analogue à (fTUl) . Mais l'expression de Vx(a) est telle qu'elle ne permet pas 
d'estimer la fonction / à toutes les fréquences (sauf éventuellement en et en 00 par un choix asymptotique 
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4.2 Analyse par ondelette et estimation de la densité spectrale à partir d'une trajec- 
toire observée à des instants discrets aléatoires 

Pour une trajectoire obtenue à des instants discrets (ifc)o<fe<n) on définit un estimateur non-paramétrique 
de la densité spectrale en remplaçant les coefficients d'ondelette théoriques par les coefficients d'ondelette 
"discrétisés" selon la formule (lÏÏl). On en déduit 



Les TLC et ([T2]) restent vrais, mais les vitesses de convergence sont plus faibles: 




(13) 



Théorème 4.1 Soit X un processus gaussien défini par HP ayant une densité spectrale f vérifiant (Fl) 
et (F2). On suppose que ip est une ondelette vérifiant VF(1,4, 1/2) avec suppip C [—A, A], et X est observé 
selon échantillonnage vérifiant l'hypothèse (S(s)) avec 3+ (2H — ^ , jff — f)_|_ < s < cîo. De plus, s'il existe 
a G]1/3,1[ et 

. sis = oo,nôl+^ etr„n-(i+^)/(2+^^) 0, 

ra— >oo n— >oo 

. Sis<00,nôl''"~^^^ Q ^t^^n-{siH+l)+2a)/(s(H+2)^l+a(l~H)) _^ 

n^oo n— >oo 

i) Le TLC U^) est vérifié également par /ra(Ç); 

ii) Pour toute bande de fréquences finie < luq < < oo, on a le développement limité du MISE: 



]E 



5 Applications numériques 

Résultats de simulations pour des mouvements browniens fractionnaires 

La densité spectrale d'un mouvement brownien fractionnaire de paramètres {H, a"^) s'exprime sous la forme 
fH,a^{C) = t:~^HT{2H) sm{7rH)a'^ l^]^^^"^^ . Nous avons simulé des trajectoires de différentes longueurs 
de mbf pour différentes valeurs de H, avec a'^ = tt {HT{2H) sin(7riî)) , de telle manière que la densité 
spectrale soit Il est à noter que dans une telle simulation le choix de a s'avère crucial: si a est 

choisi trop proche de la valeur frontière 1/3, l'algorithme d'estimation n'est pas stable; de même si a est 
trop proche de 1, l'estimateur ne converge pas assez vite et sa variance est très importante. De multiples 
simulations nous ont montré que plus N est grand, plus a pouvait être choisi proche de 1/3. La Figure 1 
donne un exemple d'une telle estimation pour N = 20000, II = 0.2 et a = 0.6. 



Application aux fréquences cardiaques instantanées pendant un marathon 

Les données que nous considérons ici ont été obtenues par le Laboratoire INSERM LEPHE, dirigé par la 
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-0.5 
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Figure 1: Estimation de la densité spectrale (à gauche) et de son logarithme (à droite) pour une trajectoire 
de mbf telle que H = 0.2 et N = 20000 échantillonné suivant un pas suivant une loi exponentielle de 
paramètre 1. 

des durées entre deux battements de coeur de 50 marathoniens de bon niveau courant le Marathon de 
Paris 2004. Ces données ont également été étudiées dans Kammoun et al. (2007), en considérant la 
suite des fréquences cardiaques instantanées comme une série chronologique classique (à pas régulier). 
Ici, nous considérons l'évolution de la fréquence cardiaque instantanée en fonction du temps réellement 
écoulé (suivant un pas de temps irrégulier correspondant aux durées entre deux battements successifs). La 
Figure 2 présente un exemple sur un athlète (ayant donc couru ce marathon en 2h45mnl2s) : 




1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

Temps en secondes Temps en secondes 



Figure 2: Fréquences cardiaques instantanées (à gauche) et cumulées (à droite) d'un marathonien en 
fonction du temps de course 

Nous avons appliqué l'estimateur de la densité spectrale /n(0 aux des fréquences cardiaques instantanées 
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accroissements stationnaires considérées ici). Cependant, comme cela avait déjà été proposé dans Kammoun 
et al. (2007), on a préféré découper cette série en trois sous-séries (obtenues à partir d'un algorithme de 
détection de ruptures en moyenne et variance) représentant le début (4802 données), le milieu (19708 
données) et la fin de course (2590 données) afin de tester une éventuelle évolution de l'indice d'auto- 
similarité du signal fréquence cardiaque en fonction de la fatigue. Notons que le pas d'échantillonnage ici 
suit approximativement une loi gaussienne (donc S{s) est vérifié pour tout s > 0) de moyenne 0.376s et 
d'écart-type 0.01s. Pour appliquer les résultats précédents, on considère la base de temps en minutes et 
ainsi (5„ — 0.006mn. La Figure 3 ci-dessous représente le logarithme de la densité spectrale en fonction du 
logarithme des fréquences pour la partie de fin de course. 




Logarithme des fréquences choisies (Hz) 

Figure 3: Logarithme de la densité spectrale en fonction du logarithme des fréquences (fin de course) 

On retrouve dans la Figure 3 la même propriété évoquée dans Kammoun et al. (2007): il existe une 
phase de linéarité dans ce graphe. Ceci avait permis la modélisation de ces données par un mouvement 
brownien localement fractionnaire, qui est une généralisation du MBF (même densité spectrale dans une 
bande de fréquence) et le coefficient directeur de la droite fournit une estimation du paramètre de fractalité 
locale H. On trouve, de la même manière: 

~ 1.12 (début de course), Hm — 1-15 (milieu de course) et Hd ~ 1.38 (fin de course). 

On constate une augmentation de la valeur de ce paramètre en fin de course, ce qui donne une intéressante 
caractérisation de l'apparition d'une fatigue du coeur (voir Kammoun et al, 2007). Ce résultat devra 
être étayé par des études ultérieures afin par exemple d'expliquer l'augmentation depuis une décennie de 
l'occurence de mort subite au cours des marathons. 

Remerciements : Nous tenons à remercier Mario Wschebor (Université de Montevideo) et Antoine 
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